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Abstract 

We study a system containing many particles of identical kinematics with a zero 
range interaction that scatters one from the other, and with the possible exchange 
of an attribute. Taking an initial condition in which the attribute is asymmetrically 
distributed in the regions of momentum space occupied by the particles, we study 
the rate at which it becomes uniformly distributed, through collisions. We find, in 
some circumstances, a rate that is much faster than that which would be estimated 
from cross-sections. This behavior is attributable in some general sense to N-particle 
entanglement. We suggest applications to neutrino physics, where the attribute is 
neutrino flavor. 

FERMIL AB-Pub-03 /062- A 

1 Introduction 



We consider a system consisting of stable or quite long-lived particles (i.e. 
stable for the time range that in which we shall be interested), confined in 
a box for simplicity, and interacting occasionally with one another. Taking 
an initial state which is, at least with respect to some of its attributes, not 
in statistical equilibrium, we can discuss time scales for evolution of the gross 
features of the system. In general these scales are determined by cross-sections. 
But an exception to this assertion about scales and cross-sections can be found, 
among other places, in the transformation of one state of a particle into another 
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through evolution that can occur in isolation from other particles, such as in 
neutrino oscillations or the precession of spins in external magnetic fields. In 
this case it may not be to the point to discuss the instantaneous "rates" at 
which the system changes, since we are watching coherent development of an 
attribute of the system with time behavior generically like sin cut. 

In the present note we demonstrate some circumstances under which we do 
not have such slow single-particle "precessions" but in which there is a long 
term coherent process arising from the short range interactions of particles, 
which can be randomly distributed in momentum space. This behavior can 
arise when the particles have another attribute that can be traded in the 
course of an interaction. Both to have a concrete framework and with a view 
to a possible application, we shall consider a system of neutrinos. In order 
to eliminate extraneous effects, we consider an example in which there are 
two flavors of massless neutrinos, designated i/e and 1/^, and there is no neu- 
trino flavor mixing in the Hamiltonian. There is an initial distribution of the 
neutrinos specified, one in which there is some systematic fiavor asymmetry, 
such as neutrinos of one flavor being predominantly of higher energy than 
those of the other flavor, or with different angular distributions for the two 
species. The question we pose is: "At what later time would the distributions 
would become more or less equal through neutral current interactions?" The 
conventional answer, as mentioned above, is that this time is determined by 
a scattering rate proportional to the weak cross-section times the density of 
scatterers. Equations in the literature do indeed predict shorter time scales for 
effects arising from neutrino-neutrino interactions in the case in which there 
is flavor oscillation built into the Hamiltonian [1] . The effects studied in these 
papers are analogous to the "forward scattering", or "index of refraction" 
terms familiar in the study of the passage of neutrinos through matter that 
contains electrons, in that they have an inverse time scale that is proportional 
to GfP where p is the neutrino number-density of the medium and Gp is the 
Fermi constant. This is in contrast to cross-section effects, which are of order 
G'j^^puj'^, where uj is the energy of the scattering particles. These studies, how- 
ever, predict that such effects are strictly absent, for pure neutrino systems, in 
the absence of flavor mixing in the Hamiltonian, if one one begins, as we do, 
with a state that is diagonal in the flavor space [2]. The essential difference 
between our approach and that of previous work is that we retain much of the 
full complexity of the multi-body physics involved, rather than assuming we 
may describe our ensemble of particles with a single body density matrix. ^ 



Some of the subtleties involved in neutrino-neutrino forward scattering have re- 
cently been re-examined by A. Priedland and C. Lunardini [3]. In particular, they 

raise the point that there may be circumstances in which a single body description 
is inadequate, as is certainly the case for the effects we study here. However, we 
do not see how to obtain the results of the present paper using the perturbative 
approach of these authors. 
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2 Model 



The basic process that we consider is simply 

i/e(p) + i^r(q) i^e(q) + i^r(p), (l) 

where the momenta p and q are drawn from the initial distributions. These 
processes will tend to reduce the correlations of flavor and momentum that 
we assume are present in the initial state. The question is: "At what rate?" 
The part of the Hamiltonian that will provide our effects operates only in the 
subspace of the initial momentum states. We define annihilation operators 
for a i^e of momentum and bi for a i>t of momentum pj, where i runs from 
unity to the number of single particle momenta. ^ Given an initial state of 
A^i momenta occupied by i/g and N2 momenta occupied by i/t, we take the 
effective Hamiltonian that implements the full set of reactions in Eq. (1) to 
be, 

Hi'"' = E h + a|a,6]6.] (2) 

where V is the volume of the system and the weight function fij is of or- 
der unity. The sum extends over all of the {Ni + N2) momentum states 
of the system that are initially occupied by either flavor of neutrino. We 
have omitted the terms in the Hamiltonian corresponding to the processes, 
^e{Pi) + ^riPj) ^e{Pi) + ^APj), which do uot Contribute to our effects. In- 
cluding such terms, that is to say, a diagonal contribution to the Hamiltonian, 
alters the wavefunction of the system only by an irrelevant overall phase. Since 
the energies of the basic set of unperturbed states are exactly the same, the 
time evolution of the system will be entirely determined by this interaction 
Hamiltonian. 

In a realistic problem the form of fij will depend on the circumstances of 
the application. For example, if we had started with the complete form of 
the neutral current, neutrino-neutrino interactions, then the matrix elements 
of the Dirac matrices in the V,A structure dictate a factor of [1 — cos(^p.^p^.)], 
where the angle is that between the two momentum vectors labeling the states. 
For simplicity, we have taken fij = 1 for some of our analytic estimates. 
In the numerical calculations we have taken a distribution of values in the 



^ Not necessarily equal to the number of particles, since in the initial state a v^, and 
a i^r could both have the same momentum. 
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range 0.5 — 1, given by 



fij = 0.5 + 



0.5 



(3) 



3 N+N 

In our first example we shall take a set of 2N momentum states to be occupied 
half by UgS and half by i/^-'s. For the initial state we take the first (bottom) A'^ 
states to be filled by z/^'s and the last (top) N states to be filled by UgS. We 
refer to this state as |^o)- Explicitly, 

|*o) = We{Pl) ■ ■ ■ T^e{Pi) ■ ■ ■ MPn)MPN+i) ■ ■ ■ J^APj) ■ ■ ■ T^t{P2n)) ■ (4) 

The interactions given in Eq.(l) exchange energy and momentum between 
particles of the two different species. An example of such a process is the 
interchange of our initial state with 

ke(Pl) ■ ■ ■ ^riPi) ■ ■ ■ l^e{PN)MPN+l) ■ ■ ■ l^e{Pj) ■ ■ ■ MP2n)) ■ (5) 

The total set of states of the complete system that we have to deal with, |^a)) 
are the n,, = (2A^)!/(A^!)^ distinct states in which the flavor indices in the 
initial state are permuted within the defined subset indexed by pi. We shall 
adopt an ordering of these states such that in each of the first ns/2 states in 
the list, the top state in the single particle list is occupied by Ue- Having taken 
all of the N of the z/g's on the top, in the initial state, we wish to estimate the 
time for "equilibration" , in the sense of the 's being more or less distributed 
with 50% in the bottom set of states. 

In the Ug dimensional space spanned by our basis states, |^a); the effective 
Hamiltonian matrix as determined from Eq.(2) has iV^ off-diagonal elements 
in each row (or column); we denote this matrix by M . Each of the UgN'^ 
off-diagonal elements of M is one of the 2N(2N-l)/2 values of the function fij 



for each one of the basis states |^q). The square root of the average squared 



of Eq. (2). 



Taking fij = 1, we have. 



(6) 
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energy of the Ug eigenstates is thus given by 



1/2 



-\/2G 



FP 



(7) 



where the trace operates in our space of Us states, and we have substituted 
the total (z/g + Vr) number-density, p, for 2N/V . We denote the corresponding 
set of eigenstates of M by I^'e), where the E^s stand for the eigenvalues. 
The state of the system at time T, which we denote by \^{T)) is given by 



|*(r)) = £|^'^;)(*s|*o)e- 



-iET 



(8) 



From this we can compute, for example, the probability P (persistence) that 
a particular one of the top 715/2 states, which was occupied at time T=0 by a 
Ve-i is occupied at time T by a Ur- For example, with the |\l'o) and the ordering 
of the Us states as described above, we ask for the probability of continuing 
to have a in the very top state, obtaining 

ns/2 

i^=El(*a|*(r))r. (9) 

a=l 



where the sum is over all the ns/2 states that have a on top. 

Of course at the moment we have said nothing about the coefficients {^a\^ e) 
that enter in Eq. (9) and Eq. (8). We anticipate that an eigenvector of a 
matrix such as M typically projects significantly onto many of the basis states 
\^a)- Our numerical tests appear to sustain this conclusion, which we will not 
quantify further in this note. Then we conjecture that the time in which P 
becomes significantly less than unity on the average, i.e. the mixing time, 
is determined by the magnitude of a typical eigenvalue. From Eq. (7) we 
already found that the root mean square of an eigenvalue is of order pGp. If 
a macroscopic fraction of the eigenvalues, that is, a number A Ug of them, are 
of this order, ^ then the effective mixing time should be of order in 
view of the relations Eq.(8) and Eq.(9). This is the same as the typical time 
scale for what the "index of refraction" or "forward-scattering" effects would 
be in the case of in which we replaced the neutrino density by an electron 
density of the same magnitude, a time scale much shorter than any effect of 
nonforward scatterings (which scale as G'p^.) 

The distribution of the eigenvalues can be determined analytically in the case 
in which fij = 1. Consider a system in which the number of Vf. and Vr is given 

^ Of course, if essentially all of the strength were concentrated in a number of 
eigenvalues which did not grow as n^, then we would not obtain macroscopic effects. 
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by Ni and N2 respectively. We shall assume, without loss of generality, that 
^1 ^ ^2- The distribution of eigenvalues have a pattern, from which we can 
observe that there are A^^i + 1 distinct eigenvalues 



V2G1 
V 



■ [{N,-t){N2-i)-i], 



0,1... iVi, 



(10) 



with degeneracies given by 



i\{Ni + N2 



(ATi + ATs)! 



l)\{Ni + N2 



1))! 



:iii 



In Fig. 1 we plot the distribution of these eigenvalues for the present case of 
N each of i^g and i^t, in the limit of large N. On the same figure we show the 
distributions obtained numerically, in the case that fij is given by Eq. (3), for 
the cases 2A'" = 12 and 2A'" = 14. Since the matrix M is traceless, the eigen- 
values sum to zero. One can see from Fig. (1) that the bulk of the eigenvalues 
(say, > 95%) he in the range 



1 Ei V 

- < — ^ — 

2 - 2A^ y/2GF 



<+l, 



;i2) 



and thus the typical energy difference between a pair of eigenstates will be of 
order (^/2Gf/V){2N) = ^/2Gpp. The qualitative features of the eigenvalue 
distribution, are not affected by taking a distribution of fi/s. 




Fig. 1. The distribution of energy eigenvalues for a system of N v^s and A'' Vt^s, 
where the energy is in units of \/2GfP- A point on the curve represents the fraction 
of states with energy below E. The heavy curve is the Umiting result for large 2N 
for the case fij = 1. The heavier and lighter dashed curves are the results for 2N=14 
and 2N=12 respectively, with fij given by Eq. (3), but multiplied by a factor of 4/3 
so that the average coupling strengths in the comparisons are the same. 

We also make a limited computational check on our heuristic estimates of 
the effects on rates, by directly solving the Schrodinger equation, beginning 
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with the above initial state, Eq. (4), and determining the probabihty that a 
particular one of the z/g's in one of the upper states has been replaced at a 
later time by a Vt. In the case 2N = 14 we have to solve = 14! / (7!)^ = 3432 
coupled linear differential equations, which is our computational limit. The 
results of these calculations are shown in Fig. 2 for a case in which the couphng 
function fij is given by Eq. (3). We see that the curves show very similar 
behavior, and except for the smallest case with 6 particles, the times elapsed 
to the first minimum are very close to each other. The results appear to sustain 
the analysis, given above, for a mixing rate that is independent of particle 
number for fixed density. Note, however, that the locations for, say the first 
peaks in the respective curves move to greater time as N is increased, albeit 
at a slower and slower rate. One cannot absolutely conclude from these data 
that there is a limiting point for large N and fixed density. Provisionally, 
taking the mixing time to be that in which the curves cross in Fig. 2, we 
find Tmix ~ {V^Gfp)~^ in the limit of an infinite system. Unfortunately, a 
perturbation theoretic approach is valid only for small times (smaller and 
smaller times as the particle number is increased) and cannot shed any light 
on the outcome. 

1 

0.8 
0.6 

P 

0.4 
0.2 



5 10 15 20 25 30 35 

T 

Fig. 2. The mixing parameter, or single level persistence probability, P, as defined 

in text, as a function of time, for a system of N and N v-t-. The curves correspond 
to the number of particles: 2A^=6, 8, 10, 12 14, where the heaviest curves are for 
highest N. The unit of time is {V^Gp p)"^ 



4 1+(2N-1) 

Next we look at a problem that is so closely related to the above one that 
one might (erroneously) guess that the behavior is similar. We start with 
the same Hamiltonian Eq. (2), and set of 2N participating states but now 
take an initial condition with only the top state filled with a Vf, and all the 
remaining (2iV— 1) single-particle states filled with z/^'s. The effective subspace 
of system states that are connected together is now Ug — 2N dimensional, and 
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is described by the location of the single z/g. But now, in the analogue of 
Eq. (6), for the case of fij = 1 we find that the average squared eigenvalue is 
(2iV)-iTrM2 = {2N-1){^/2GF/Vf, in contrast to the value, N^, obtained in 
the {N, N) case. Furthermore, there is a single state, with eigenvalue (2A'" — 1), 
while the remaining (2A'" — 1) states all have eigenvalue —1. Denoting the state 
with large eigenvalue by ^5, we note that | (\E'5|\E'o) P = so that the effects 
on P of mixing with this state are of order N^'^; and since there are no energy 
differences among the remaining A^— 1 eigenstates, the total effect (at the order 
Gf level) will vanish in the limit of large N. We have confirmed numerically 
that the introduction of a scatter into the coupling constants, using Eq. (3), 
does not change these conclusions qualitatively. The results for persistence 
versus time in this case are shown in Fig. 3, for the case of several values of 
A^. They clearly show the behavior in the short and intermediate time 
regions. 

1 

0.9 
0.8 

P 

0.7 
0.6 



20 40 60 80 100 

T 

Fig. 3. The same as Fig.2 except for the case of an initial system state with only 
one particle state occupied by a z/g and the others occupied by z^r's, with P the 
probability that this the is in its initial state. We compare the cases 2A'^=30, 60, 
120, 240 where the slower evolution corresponds to the higher number of particles. 



5 Discussion 

Thus far in our discussion, we have used "z/g" and "z/^" merely to label our 
states, and we have given no hint as to how the above behaviors could make a 
difference to observable results in the neutrino world. In fact there are two sit- 
uations that have been discussed in the literature in which neutrino-neutrino 
scattering is thought to play a role, and in which the size of the parameters 
is exactly such as to make effects with our time scale {\/2Gfp)~^ of possible 
importance, and indeed as important as the effects noted in the literature. The 
first is in early universe scenarios in which one has assumed some degree of 
neutrino degeneracy [4]. The second is in the region just above the neutrino- 
sphere in the supernova process, where differences in the energy spectra for 
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different neutrino flavors play a key role in whether or not there is an efficient 
R process for nucleosynthesis [5]. Another example where neutrino-neutrino 
forward scattering is relevant, is the flavor conversion of high energy neutrinos 
from astrophysical sources as they propagate through the relic neutrino back- 
ground [6] . We emphasize again that up to this point we have not included the 
neutrino oscillation physics that drives the results in these works. Wc intend 
to address the combination of the two themes in a subsequent publication. But 
we can comment on what we believe to be the insufficiency of the previous 
works on the subject. What we have computed above, the probability that a 
specific state, say Pi, that was originally occupied by a 1/^ (in the initial state 
for the complete system I^E'o)) is still occupied by a i/e at time t, has a simple 
expression in terms of the Heisenberg picture operator, ai{t), 

P^{^o\alit)aiit)\'^o), (13) 

where the Heisenberg operators are chosen to coincide with the Schrodinger 
operators that were introduced in Eq. (2) at t = 0. Now the operator in 
Eq. (13) is a neutrino density operator for the particular state Pi. If we write 
an equation of motion for this operator, by taking the commutator with the 
Hamiltonian Eq. (2), the right hand side is a quartic in the operators a\ a, b\ b, 
and there is a sum over one index j. When we encounter, say, a term like 

^al{t)b,{t)bl{t)a,{t) (14) 
j 

we would like to be able to replace I]j(aj6j) by a density matrix element, 
off-diagonal in the fiavor space, which represents the entire average state of 
the medium. Then the multiplying operator fejoj would be a corresponding 
operator but for the single mode i. Indeed, this is exactly the assumption that 
yields the non- linear terms in the equations for the density matrices derived 
in Ref. [2] and used in Ref. [4,5]. 

As noted above, with our assumed flavor- diagonal initial conditions, and in 
the absence of neutrino-mixing in the Hamiltonian, such terms do not create 
any effect. Our effects are exactly due to the fact that the replacement of 
the four-operator product by a product of two expectation values, as sketched 
above, is not justified. Even if we were to assume that a kind of factorized 
ansatz would pick up the leading terms, which we believe is highly unlikely, 
we would want add the results of pairing, say, a]aj in Eq. (14). That is to 
say, we would be driven to consider density matrices that are off-diagonal in 
momentum space, as well as in fiavor space. 

To summarize, we have evidence for a new, macroscopic quantum effect that 
could change the outcome of calculations in which i>-i> scattering matters. The 
caveats that must be added are: 
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(i) Accepting the earlier conclusions for the model defined by the Hamiltonian 
Eq. (2), we must go back and ask whether we have defined the correct problem. 
All those modes that we left out, when we truncated our Hamiltonian to 
momentum states that were initially occupied by one species or the other, 
can they make a difference? These additional modes would be populated via 
non-forward scattering interactions, which enter the problem at the level of 
cross-section (x density), that is, on time scales of order {G'^puj'^)^^ . Our 
expectation is that they will not affect our results on the much shorter time 
scales we consider (of order {Gfp)~^- 

(ii) Our calculations pertain directly to plane waves in a volume V . Ordinary 
neutrino oscillation theory has come under fire repeatedly from authors that 
suspect that preparation-of-state considerations bring the idealized, plane- 
wave picture into question. These criticisms have been successfully answered 
more than once [8] for the case of single-particle oscillation, but they could 
surface again in our present context. 

At the present time, we are cautious about claiming that our results will 
be important to the neutrino physics either in the early universe or in the 
supernova. These are, however, the two most obvious applications in which 
the neutrino number density is high enough for neutrino-neutrino scattering 
to be important. We further note that the flavor-energy correlation is critical 
to understanding the physics just outside the supernova neutrinosphere. Of 
course, we would need to include neutrino mass and mixing in the model 
in order to address a realistic situation. Superficially, we can note that the 
energies attributable to mass effects are small compared to the inverse time 
associated with our processes under the conditions that prevail just outside the 
neutrinosphere of a supernova. But more analysis is required before reaching 
a conclusion in that context. A negative feature of our work is that in contrast 
to being able to use the quite simple equations for the density matrix posited 
by Ref. [2], we appear to be doomed to treating the full complexity of the 
many-body physics that arises in these systems. 



6 An example with a finite number of discrete states 

We demystify the physics of much of the above, to some degree, by thinking of 
a fairly large but finite set of discrete quantum subsystems with interactions 
of comparable strength between each pair. In addition to being an interesting 
exercise this also is relevant to applications in other arenas than the world of 
neutrinos. We consider a set of N spins associated with one set of fixed site 
labels, say red, and another set of N spins associated with a second set of 
labels, say blue. The interaction exchanges all pairs of spins with comparable 
strength, or order for each pair. Let us now start with all reds up and all 
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blues down, and consider the short time evolution of this state. The probability 
that the spin on a particular red site becomes down in a very short time t is 
of order h = g^Nt^. We might think that this perturbation calculation stays 
more or less valid until such a time that h is, say, 0.1. But if we go back 
through the same mechanics that we used for our continuum problem, we see 
that is not be the case; the system will mix in a time T such that gNT = 1 ^ , 
at which time the above perturbative estimate of the amount of mixing is only 
h = 1/N. 

Next consider the case in which the first red spin is up, and every other red 
and blue spin is down. In complete correspondence with the (1, 2N-1) neutrino 
case, we find that the above perturbative estimate of the mixing time, for the 
first red spin to become appreciably blue, is correct. The mixing time for this 
single state is longer than that of the first example by a factor of N, for large 
N. 

What makes these cases so different from each other? In both cases the spin 
at the first red site interacts with each of the N blue spins. But in the second 
example, the blue spins themselves interact only with that one red site, at 
the turn-on time, while in the first example every blue spin is being affected 
by every red spin from the beginning. The difference lies not in the number 
of interactions that the red spin sees, but in the entanglements of the states 
that it interacts with. Since in the first case we get faster evolution than we 
would have expected from the perturbation estimate we classify the effect as 
a "speed-up" . And the reason is clearly the multiparticle entanglement within 
the system. For completeness, we should note that even in the second example, 
the blue sites, with which the spin on the distinguished red site interacts, do 
develop mutual entanglement, but with a much slower initial rate than in the 
first example. In effect, they see each other only through their mutual coupling 
to the distinguished red site. 



7 Conclusion 

We have examined a "speed-up" of evolution through entanglement, both in 
a discrete system of spins and in a neutrino model with a continuum of states. 
Though the practical meaning of "speed-up" is slightly different in the two 
cases, the formal source of the effect is the same. 

The reader might ask the question: "Are the authors using the word 'entan- 

^ Note that in this case we do not redefine the couphng constant with increasing 
N; there is no volume or density involved in the models. Thus the eigenvalues grow 
as N. 
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glement' in some precise mathematical sense, or are they just using it to say 
the system gets comphcated?" We contend that it is the former, in the sense 
that the state of the system cannot, in general, be factorized into a product 
of single particle states. Note, however, that the quantification of multipar- 
ticle entanglement is not a concept that that has been precisely defined in 
the literature. Definitions of two-state entanglement, however, have received 
much attention. The reader can consult Ref. [7] for a demonstration of how 
in the simplest relevant ordinary Schrodinger example, with two particles in 
a double well, entanglement in a precise mathematical sense is spontaneously 
generated if the two particles interact with each other. 

Returning to the neutrino example, the first order effects that mixed the states 
had an inverse time-scale ^j2Gpp\ this is the "speeded-up" rate. The usual 
inverse time-scale for (non-forward) scattering effects is of order G^piJ^ where 
Lo is of the order of the particle energies and much, much slower, than the 
"speeded-up" rate when lo is of order an MeV. Comparing again with the 
calculation with only one i^e in a sea of i^t-'s, where there is no speed-up, the 
difference is attributable to the degree of entanglement. 

The inclusion of entanglement requires abandoning the single-particle descrip- 
tion of the system, to include the full complexity of the many-body physics 
involved. Single-body descriptions do predict significant neutrino-neutrino for- 
ward scattering effects [1,2], which, however, are absent for the flavor diagonal 
initial states we consider here. The results of our many-body calculation may 
thus have interesting consequences in situations where v-v scattering is im- 
portant. 
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NOTE ADDED 

After our paper appeared, Friedland and Lunardini [9] analyzed the totally 
symmetric case, where all the coefficients jij = 1, and have shown that the 
mixing times (for couplings scaled to 1/ [particle number]) grow as \/N, for 
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large N, in contradiction to our conjecture. However, this analytic solution 
cannot be readily extended to the case of unequal fij, so the results of [9] are 
perhaps inconclusive in this case. As explained in the text, the simulations 
presented in Fig. 2 are for the case of scattered coupling strengths. These 
simulations are not adequate for distinguishing \/N limiting behavior from 
constant limiting behavior, due to the restrictions on N imposed by computa- 
tional resources. 

We have since applied the techniques of [9] to some cases of differing coupling 
schemes and initial conditions, but still with sufficient symmetry to do simu- 
lations for much larger values of N than previously, and we find evidence for 
fast evolution in some of these more generalised systems. 
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